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Abstract. A simple method for describing the evolution of a quantum state of a double 
quantum dot system interacting simultaneously with the electromagnetic environment and with 
the lattice modes is developed. It is shown that the joint action of the two reservoirs leads to 
nontrivial effects in the system dephasing. As an example, the impact of phonon-induced initial 
dephasing on the radiative decay of delocalized exciton states is discussed. 



1. Introduction 

Double quantum dots (DQDs) are systems composed of two quantum dots placed close enough 
to each other for new coherent and collective phenomena to occur which cannot be reduced to 
the known properties of individual dots. Such new phenomena related to the optical properties 
of the system include modification of the optical response due to inter-dot coupling [if; 0] or 
collective effects in the spontaneous emission 0; 0]. Another class of effects is related with 
carrier-phonon interaction. Examples include phonon-assisted exciton transfer between the dots 
@; S 0]' modification of optical spectra due to phonon packets traveling between the dots 
or phonon-induced decay of entanglement 0|. The special properties of DQDs may open the 
way to new applications, like long-time storage of quantum information [lOj], conditional optical 
control [l[ that may lead to an implementation of a two-qubit quantum gat e Hill, g eneration of 
entangled photons [12] or coherent optical spin control and entangling [13l ; Il4l: 15]. Although 



both spontaneous emission and phonon-related phenomena have been studied to some extent 
for DQD systems, the interplay between these two classes of effects has not been investigated. 

In the present contribution, we formulate a theoretical description of the evolution of a double 
quantum dot coupled simultaneously to the two reservoirs: the quantum electromagnetic field 
(photon vacuum) and lattice vibrational modes (phonons). After defining the model (Sec. |2J), 
the approach to the simulation of the quantum open system is presented (Sec. [3]), followed by 
an example of an interplay between the two couplings (Sec. SJ). 



2. The system 

The system under study is composed of two stacked self-assembled QDs with transition energies 
ei and 62, interacting with their phonon and photon (radiative) environments. We restrict the 
discussion to the ground states of excitons in each dot and assume that the spin polarizations 
of the excitons are fixed. As the exciton dissociation energy in absence of external electric fields 
is rather large (several to a few tens of meV), we consider only spatially direct exciton states, 
i.e., such that the electron-hole pairs reside in one and the same dot. 



Thus, the model includes four basis states: the state without excitons, |0), the states |1) 
and 1 2) with an exciton in the first and second dot, respectively, and the 'molecular biexciton' 
state 1 3) with excitons in both dots. The system evolution will be described in a 'rotating basis' 
defined by the unitary transformation 

U _ e iEt(|lXl|+|2)<2[+2|3X3|)/ft 

where E = (e\ + 62)/2- Note that this is not equivalent to interaction picture with respect 
to the Hamiltonian of uncoupled dots. Upon the standard weak-coupling derivation of the 
Master equation, the latter would always lead to uncorrelated emission for non-identical dots. 
On the contrary, the present transformation allows us to explicitly keep the energy difference 
A = (ei — €2)/2 which yields unitary evolution superposed on the dissipative one, consistent 
with the result of the Weisskopf-Wigner approach [3J], and allows one to correctly account for 
the transition from "identical dots" to "different dots". The Hamiltonian is then 



H — -f^DQD + H ph + H vaA + H c _ ph + -f^c-rad- 

The first term describes exciton states in the DQD structure, 

H DQD = A(|1}(1| - |2)(2|) + V(\l)(2\ + |2)(1|). 



(1) 



where V is the coupling between the dots, which may originate either from the Coulomb (Forster) 
interaction or from tunnel coupling (it can be assumed real). It is convenient to introduce 
the parametrization A = £cos2#, V = £sin2#, where 28 is the energy splitting between 
the single-exciton eigenstates of -£Tdqd and 9 is the mixing angle of the single-exciton states. 
Electron and hole wave functions are modelled by identical anisotropic Gaussians with identical 
extensions I in the xy plane and l z along z for both particles, 



V>i,2 (r) ~ exp 
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where D is the distance between the dots. 

The phonon modes are described by the free phonon Hamiltonian H p h = ^ fc hu>kb\,bk, where 
bu are bosonic operators of the phonon modes and u>k are the corresponding frequencies. 

Interaction of carriers confined in the DQD with phonons is modelled by the independent boson 

Hamiltonian 
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where fi ' are system-reservoir coupling constants. For Gaussian wave functions, the coupling 
constants for the deformation potential coupling between confined charges and longitudinal 
phonon modes have the form f^' 2 ^ = /fce ± * fcz ' D / 2 , where 



fk 
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Here v is the normalization volume, k±/ z are momentum components in the xy plane and 
along the z axis, a e /^ are deformation potential constants for electrons/holes, q is the speed of 
longitudinal sound, and g is the crystal density. We assume that off-diagonal carrier-phonon 
couplings are negligible due to small overlap of the wave functions confined in different dots. 



The third component in our modeling is the radiative reservoir (modes of the electromagnetic 
field), described by the Hamiltonian -£f ra d = x^ w k c k \ c k,x, where ci^Xi c h\ are photon 
creation and annihilation operators and are the corresponding frequencies (A denotes 
polarizations). The QDs are separated by a distance much smaller than the relevant photon 
wavelength A = 2irhc/E, so that the spatial dependence of the EM field may be neglected (the 
Dicke limit). The Hamiltonian describing the interaction of carriers with the EM modes in the 
dipole and rotating wave approximations is 

tf c _ rad = £_ Yl e ~ iEt ' h 9kxA,x + H - c - . ( 3 ) 

k,X 



withS_ = |0)(1|+|2)(3|+|0)(2|+|1)(3| and g kX = id-e x {k) J ^g^, where d is the interband dipole 
moment, £q is the vacuum permittivity, e r is the dielectric constant of the semiconductor and 
e\(k) is the unit polarization vector of the photon mode with the wave vector k and polarization 
A. For wide-gap semiconductors with E ~ 1 eV, zero-temperature approximation may be used 
for the radiation reservoir at any reasonable temperature. 



3. The system evolution 

In certain limiting cases, analytical formulas for the evolution of the DQD system may be 
found. For uncoupled dots (V = 0) interacting only with lattice modes (phonons), an exact 
solution is available 0. If only the radiative decay is included, a solution in the Markov limit 
can be obtained 0]. Here, we propose a general description which allows one to deal with 
the simultaneous action of both these environments. We describe the evolution of the reduced 
density matrix of the DQD system in the interaction picture with respect to i^DQD by the 
equation 

P = £radb] +£ph\fi\- 

Here the first term describes the effect of the radiative decoherence in the Markovian limit in 
terms of the Lindblad dissipator 



£rad[/?] = F ra d 



S_(t) /0 S + (t)-i{S + (t)S_(t),p} + 



where E_(t) = Ej_(t) = e <HDt »*/ R E_e -ifl ' Dt » Dt / ft and Frad _ E 3 \d\ 2 ^T r /(3Tre c 3 h A ) is the 
spontaneous decay rate for a single dot. The second term accounts for the interaction with 
the non-Mar kovian phonon reservoir. We use the time-convolutionless equation 



£ph[p] = ~ drtr ph [H c _ ph (t), [H c - ph (r), p(t) ® p ph \] , (4) 
J o 

where i? c _ p h(t) = e t ^ Huc ^ +H ^ t / h H c ^ p \ 1 e^ H ^ c ^ +H ^ t / h is the carrier-phonon interaction 
Hamiltonian in the interaction picture, p-^ is the phonon density matrix at the thermal 
equilibrium, and tr p h denotes partial trace with respect to phonon degrees of freedom. While 
some phonon-assisted processes allow a Markovian limit, there are many effects that can only be 
reproduced if the reservoir memory is included (such as the initial phonon-assisted dephasing) . 
Moreover, there seems to be no universal way to extract the Markov limit in various physical 
situations. By using the general Eq. @, we obtain an all-purpose equation of motion at a 
modest computational cost. 

In the evolution generated by the Hdqd, the states |0) and |3) are invariant and nontrivial 
evolution takes place only in the subspace spanned by |1),|2). The corresponding evolution 



operator can therefore be found easily. Transforming Eq. ([2]) to the interaction picture and 
substituting to Eq. we find the phonon-related contribution to the system evolution in the 
explicit form 

VM = E [ T i(t)p(-t)Si(t) - Si(t)Ti(t)p(jt)]+n.c. 

i=l,2 

Here the nonzero elements of the matrices Si are 

1 / 2ft\ 

<l|5i|l) = 1 - <2|5i|2) = 1 - <1|5 2 |1> = <2|5 2 |2) = -sin 2 20 1 - cos — , 



h J 

(l|5i|2) = <2|5i|l)* = -<1|5 2 |2> = -(2\S 2 \1)* = - sin20sin — - - sin 40 ( 1 - cos " M 
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<3|5i|3) = <3|5 2 |3) = 1, 



and the matrices Ti(t) are defined by 

T i(*) = Y, f drS^Rijit 
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Here 
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Rn(t) =R 2 i(t) = ^^l/fcpcoB^Dfnfce^ + Cnfc + ^e-^**] 

are memory functions with nfc denoting the Bose distribution for phonon modes. The above 
equation of motion for the reduced density matrix strictly reproduces the results in the limiting 
cases mentioned above. Moreover, for the case of a DQD coupled to phonons with non-vanishing 
inter-dot coupling V, it yields results reasonably close to those obtained by a much more complex 
correlation expansion technique jjj. 

In numerical simulations, we take the parameters corresponding to a self-assembled 
InAs/GaAs system: a e — = 9 eV, p = 5350 kg/m 3 , q = 5150 m/s, the wave function 
parameters I = 4.5 nm, l z = 1 nm, D = 6 nm, and the radiative recombination time (for a single 
dot) l/r rad = 400 ps. 



4. Example: Impact of phonon-induced pure dephasing on spontaneous emission 

As an example of an application of the formalism developed above, let us study a simple 
example of interplay between phonon-related and radiative phenomena. First, we consider 
identical, uncoupled dots, that is, V = 0, A = 0. In such case of two identical emitters, the 
spontaneous emission has a collective character which leads to a strong modification of the 
exciton recombination process 0]. In particular, out of the two delocalized single exciton states 
|^±) = i |2))/V% \ip+) has a "superradiant" character and decays twice faster than a 
single-dot excitation, while is stable against recombination ( "subradiant" ) . In Fig. [TJa), 
the decay of the state in the absence of phonon perturbation is shown by the black dotted 
line. 

For a single dot, phonon-induced dephasing (within an independent boson model with a 
harmonic reservoir) leads only to coherence decay on picosecond time scales and does not affect 
the radiative recombination on long time scales. This situation changes in quite an interesting 
manner in the two-dot case. As one can see in Fig. [D^a), in the presence of phonon-induced 
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Figure 1. Decay of the exciton occupation for the superradiant initial state in the 

absence of phonon-induced dephasing (black dotted lines) and in the presence of phonon-induced 
dephasing at various temperatures as shown, (a) identical, uncoupled dots, A = V = 0; (b) 
non-identical, coupled dots, A = V = 1 meV. 



dephasing, the exciton recombination slows down and a long-lived occupation appears in the 
system. Moreover, the final exciton occupation grows as the temperature increases. Obviously, 
this effect must be due to some kind of coherence leading to occupation trapping in a subradiant 
state which is only possible in such a clear form because the model does not contain dephasing 
channels that would destroy this coherence (like anharmonicity- induced scattering fl6l]). Still, 
such a noise-induced coherence is remarkable. 

The reason for this behavior is that during the first few picoseconds of the open system 
evolution the coherence between the single exciton states decays due to carrier-phonon dynamics 

h 



according to (l|p|2) = \e h ® , where p is the density matrix of the exciton subsystem and 
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where i p h is the characteristic time of the phonon-induced initial dephasing which is of the order 
of a few picoseconds [9J. Actually, the above formulas are valid in the absence of spontaneous 
emission. However, on the time scales t <C T, when the dephasing takes place, radiative 
contribution can be approximately neglected. As a result, the initial density matrix p = )(?/>+ 1 
is transformed into 

p' = — 2 — hMGM + — 2 — 

This dephased state is a mixture of a superradiant and a subradiant component. The latter 
does not decay radiatively and leads to the persistent occupation tail visible in Fig. [TYa). 

In the more realistic case of dots that differ in their transition energies and are coupled 
[Fig. Kb)], the radiative recombination process still depends on temperature but this dependence 
is reversed: Now, the emission speeds up as the temperature is increased. In any case, however, 
the decay in the presence of carrier-phonon coupling is slower than in the absence of phonons. 
In order to explain this behavior, we note that although the eigenstates of -ffDQD in the 
presence of energy mismatch are not purely sub- and superradiant one of them still decays 
faster than the other [3]. For V > 0, the lower eigenstate has a partly subradiant character. 
Moreover, the coupling between the dots enables excitation transfer between these eigenstates 
0| . Coupling to phonons provides a thermalization mechanism for the occupations of the two 
eigenstates. At T = this means a transfer to the lower-energy, subradiant state, which strongly 
suppresses emission. At higher temperatures, however, the contribution from the higher-energy, 
superradiant state increases and the radiative decay speeds up. 



5. Conclusion 

Double quantum dots not only show a richer structure of exciton states than a single dot but also 
admit much more complicated dephasing channels. The approach to numerical simulation of the 
quantum open system dynamics of a single DQD structure presented in this contribution opens 
the way to a systematic study of these effects. Although much more accurate (sometimes even 
exact) methods can be used in the case of two dots interacting with a single reservoir 0; 0; E3] 
treating the evolution under joint action of both reservoirs at the same level is either impossible 
or very demanding numerically. The method presented here, based on the TCL equation for 
the non-Markovian phonon-related effects and the Lindblad equation for spontaneous emission, 
yields simple equations for the complete biexciton density matrix of the DQD system and allows 
one to calculate arbitrary exciton-related properties of an undriven system (or a system excited 
with a single ultrashort pulse). 
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